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FLUIDIZED BED REGENERATORS FOR BRAYTON CYCLES 


by Lester D. Nichols 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

A recuperator consisting of two fluidized bed re- 
generators with circulating solid particles H consid- 
ered for use In a Brayton cycle. These fluidized beds 
offer the possibility of high temperature operation if 
ceramic particles are used. 

Calculations of the efficiency and size of flu- 
idized bed regenerators for typical values of operating 
parameters have been made and compared to a shell and 
tube recuperator. The calculations Indicate that the 
fluidized beds will be more compact than the shell and 
tube as well as offering a high temperature operating 
capability. 

INTRODUCTION 

The efficiency of a Brayton cycle depends among 
other things upon the turbine inlet temperature and 
how much of the heat in the turbine outlet gas stream 
can be recirculated to th ; gas from the compressor out- 
let Aj the turbine inlet temperature Increases so 
does the operating temperature of the heat exchanger 
used to recover the exhaust heat. The size of the heat 
exchanger to transfer this heat must also be kept small 
in order to keep costs low. Since fluidized beds using 
ceramic solid pal tides can operute dt very high tem- 
perature and have a large amount of heat transfer sur- 
face area in a given heat exchanger volume, they make 
likely candidates for high efficiency low cost Brayton 
cycle application. This concept was first introduced 
as air preheaters for MHD cycles. 1*2 

A schematic diagram of the Brayton cycle using 
fluidized beds Is shown in Fig. 1. The turbine exhaust 
flows into the upper of two fluidized beds with a ve- 
locity larger than the minimum fluidization velocity. 
The particles at the bottom of that bed will become 
heated. Their temperature will approach the turbine 
outlet temperature if they remain in contact with the 
gas long enough. Also, the gas will cool as it trans- 
fers its heat to the particles. The hot particles are 
removed from the bottom of the bed and cold particles 
are placed on the top of the bed to maintain a constant 
bed height. The gas leaving the fluidized bed will 
approach the temperature of the cold particles, if they 
remain in contact with the gas long enough. 

The net effect of the operatijn of this heat ex- 
changer is to transfer the heat in the outlet gas from 
the turbine to a t frzam of solid particles. The solid 
particles are flowing in a direction opposite to that 
of the gas. 

The heat in this stream of solid particles can, 
in turn, be transferred back to a gas stream in another 
fluidized bed. This second bed is the lower bed in 
Fig. 1. Again, the outlet gas from the compressor 
flows into a bed at a velocity greater than the mini- 
mum fluidization velocity. The gas is to be heated by 
the solid particles to as high a temperature as possi- 
ble, so the hottest particles are put on Cop of the 
bed. They flow down the bed, giving their heat to the 
gas until they reach the bottom of the bed where their 
temperature approaches the temperature of the compres- 
sor outlet temperature again if they remain In contact 
with the gas long enough. They are then removed from 


the lower bed and sent to the upper bed to be reheated. 

The combination of the two fluidized beds and the 
circulating solid particles constitutes the recuperator 
for the Brayton cycle. The solids from the upper bed 
could flow by gravity Into the lower bed. However, the 
lower bed is at a higher pressure (nominally the com- 
pressor outlet pressure). But the height of the hot 
particle transfer tube can be designed so as to provide 
a static "head" of solid particles equal to the pres- 
sure In the vessel Then the solids will flow into the 
lower bed. The solid particles can be returned to the 
upper bed by pneumatic transport, l.e. "dragging" the 
particles back by a high velocity stream of air. 

The efficiency of this Brayton cycle depends upon 
the amount of air required to transport the solids from 
the lower to the upper beds as well as the amount of 
air which "leaks" through the tube which Is transport- 
ing the solids from the u^per bed to the lower bed. 

The size of the two beds for a specified power 
level of the cycle depends upon the velocity of the gas 
through the beds and the height of the bed. Since the 
pressure drop of the air flowing through the beds is 
proportional to the height of the bed there Is a trade 
off between bed size and cycle efficiency. 

This paper calculates the effect of particle di- 
ameter and bed pressure drop (proportional to bed 
height) on the Blze of the recuperator and the effi- 
ciency of the cycle. Optimum values of compressor out- 
let pressure and pressure drop across the beds are 
found. 

CYCLE EFFICIENCY 
Fluidized Bed Regenerators 

The Brayton cycle which will be analyzed is shown 
in Fig. 2. While cycle performance can be Improved by 
using turbine reheat and compressor Intercooling, it Is 
not Important for the comparisons between different 
types of heat exchangers. The only difference between 
this cycle and the usual cycle is that the mass flow 
through the compressor (state points 1 to 2) is greater 
than through the turbine. The difference in mass flow 
either "blows" the particles from the lower to the 
higher, or "leaks" from the high pressure to the low 
pressure. 

The efficiency can be written as Ref. 3: 



tIGINAL PAGE IS 


i i g^T ! '■■ " ' , ii ,, i i i jp, !|?|l|^g|)j|;jpBpWWP^^ 



(Symbols are defined in Appendix A) 


and Eq. (6) becomes 


The recuperator effectiveness can be calculated 
from the heat transfer processes In the regenerators. 
Equations for gas to gas recuperators hove been formu- 
lated and solved by Hausen^ (pp. 202-213). His results 
con be used by considering (In either regenerator) the 
gas to be one fluid and the solids to be a counterflow 
fluid. 


The recuperator effectiveness can be determined from 
Eqs. (14) and (4): 


Then Hausen's equation for the lower bed becomes 


The maximum value of n , occurs when T " 1. This 
can be seen on physlcal r £rounds by noting that this 
matches the temperatures In the solids to that of the 
gas. This Is equivalent to a counterflow gas to gas 
recuperator with negligible temperature drop across 
the separating wall. In this case: 


where 


since the area of the spheres in a unit volume of bed 
Is 6(1 - eg)/ Dp (Ref. 4, p. 333). The equations for 
the upper bed are: 


The recuperator effectiveness can approach 1 
and ( (l.e. H, and H ) are large enough. 


The bed heights and the length of the hot parti- 
cle transfer line are related to the pressure drops by 
the momentum equation (Ref. 5, p. 72): 


We must evaluate the mass flow required to "blow" 
the cold particles from the lower bed to the upper bed 
The momemtum equation (Ref. 5, pp. 389-391) is: 


The solid mass flow must be equal in the two beds 
and twe transfer tubes: 


The gas mass flow is equal the same In both beds so 
that: T - T - T. 


The heat transfer coefficients are related to bed 
conditions by (Ref. 5, p.215): h D * (1/10) (PgU 0 Cp/P r ) 

Also, if there is no heat 1"et in the transfer tubes 
T 2s - T$ s and » T 5 S . Equations (5) and (8) can 

be solved for T 2s and T 3 S : 
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and that the flow area on the shell side Is equal to 
the flow area inside the tubes. 


If h|-h o *h*^ and the tube length Is L t , then 
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The heat transfer coefficient (Ref. 6, p. 547) Is: 
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It expresses the pressure required for the mass flow 
G cs . A sensitive parameter In 'his equation Is the 
void fraction, c c , in the tube. We choose the curve In 
Kunll and Levensplcl (Ref. 5, Fig. 17, p. 386). 

The "leakage" of mass thnugl the downflow line is 
(Ref. 5, p. 67): 

(p 2 - dp t + p t + ap u )PpthIP 2 ' Ap i ~ p i ' Ap u 1,lD h 
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It Is negligibly small for all cases of Interest and 
will be neglected. 

All of the terms In the efficiency equation can be 
calculated. The fluidized bed volume required can also 
be determined: 


where 
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are the terms representing the 


free board volupw. This Is the open volume above the 
bed to minimize particle carryover. We will arbltra- 
lly make the free board volume equal to the bed volume. 

The power output, P, Is: 

- (‘ * ^ @) c } <!5 ' 
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The specific volume will be dcfJr'ni a r ’» fluidized 
bed volume, V, divided by tK. yt-ve? , 

Shell and Tu b e Heat Exchanger 

The performance of the fluidized bed regenerator 
can be compared to a shell and tube heat exchanger. 
Equation (1) holds for the efficiency by setting ft 
and ft. equil to zero. The recuperator effectiveness 
is given In Jakob (Ref. 4, Eqs. 34-37 and 34-tl, 
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The length of the tubes can be related to the pressure 
drop (Ref. 6, Eqs. 21-37, p. 33): 
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The mass flow rate is: 
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where D s p Is the overall outside diameter of the re- 
cuperator. 

The power output Is obtained from Eq. (25) again 
with m), and m c - 0. The specific volume for the 
shell and tube recuperator Is: 
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SPECIFICATION OF PAR J1ETERS 


We will examine the effect: of varying certain 
parameters with respect to flu' llzed bed performance 
on the efficiency of the cycle and the size of the 
recuperators. To do this we will specify certain of 
the parameters with values typical of Brayton cycle 
systems burning a fuel and using the combustion pro- 
ducts as working fluid. 

We will choose air as a typical working fluid and 
boron nitride as a typical high temperature ceramic. 
The void fraction for the transport of the cold par- 
ticles can be found for boron nitride (density Is 
2100 kg/m3) In Ref. 5, p. 386. Typical turbine and 




coo .-•■or efficiencies and pressure drop across the 
the heat aource will be chosen. They are specified as 
follows : 


culate the specific volume of the beds by dividing the 
power output by the total recuperator volume. Uc have 
4 parameters left to vary; particle diameter, fraction- 
al pressure drop, the compressor outlet pressure, and 
the turbine inlet temperature. 


n turb 


.90 


&Pj 

P 2 


.04 


Specifying these parameters in the upward solids 
transport line equation (Eq. (22)) determines the mass 
flow rates of both the gas and the solids. This, in 
turn, can be used to calculate bed heights and recu- 
perator effectiveness and finally, to calculate effi- 
ciency. 


The Inlet conditions are T. - 330 K, 
The gas is assumed to be air and the 
evaluated at 900° K: 


properties are 


P 


4.1xl0' 5 


kg/m sec 


The system power level can be calculated by spec- 
ifying a fluidized bed dlumeter, rather than the cold 
partlclt transfet tube diameter. This will be done by 
setting the diameter of the upper (lower pressure) di- 
saster equal to its height. 


■R - 280 j/kg °K 


Shell and Tube Recuperator 


c - 1 kj/kg °K 
Pr - .7 
y - 1.4 


The same component efficiencies ond fluid proper- 
ties will be chosen as before. However, the tube di- 
ameter will replace the particle diameter. Also, the 
Mach number will be varied In addition to the pressure 
drop and compressor outlet pressure. 


Fluidized Bed Regenerators 

The solid will be taken as boron nitride with: 

p ■ 2100 kg/m 3 

c a - 1.6 kj/kg °K 

c - .973 
c 

The void fraction in the down tube Is chcsen as .4 and 
In both fluidized beds sb .45. This corresponds to 
dense and normal packing of particles with sphericity 
of .8 (Ref. 4, p. 66). The fractional pressure drops 
will be made equal |(Ap ( /p 2 ) ■ (Ap u /Pi)] and the 
height of the low pressure will be set equal to the 
diameter, H u - D u for the calculation of output 
power. If the bed height Is less than the diameter It 
may be unstable and a ''channeling" might occur through 
which most of the flow may pass. If the height gets 
much greater than the diameter then the bed may get 
unstable and "slugging" might occur (Ref. 5, p. 75). 

A parameter associated with bed performance la 
the superficial velocity (actual velocity divided by 
the void fraction). This must be higher than the min- 
imum fluidization velocity, but not so high that sig- 
nificant bubbling occurs. Higher velocity allows more 
heat to be transferred and reduces equipment size. 
However, as gas velocity Increases there will be mix- 
ing In the vertical direction because of an Increased 
number of bubbles forming and rising due to their 
bouyancy. This mixing can reduce the effectiveness of 
the heat exhehanger by bringing the Inlet and outlet 
gas temperatures closer together. Bubbles appear In 
gas-solid fluidized beds because of the large density 
ratio (Ref. 5, Eqs. (22) and (32) in Chapter 3). 
However, If the superficial velocity Is two times the 
minimum fluidization velocity (ref. 5, p. 134) then 
less than 10T of the bed volume will contain bubbles 
and mixing should be minimized. We assume that ve- 
locity ratio in our calculations. 

In order to calculate the volume of the beds we 
specify that the diameter of the pipe carrying the 
solids to the upper bed Is 10 cm. This Is important 
only in calculating the frictional pressure drop for 
the gas-tube interaction and Is unimportant compared 
to the gas-particle Interactions. Therefore, we cal- 


RESULTS 

First, a turbine Inlet temperature of 1250° K is 
selected. Then Ap/p and D„ are selected and com- 
pressor outlet pressure varied. As the pressure In- 
creases the specific volume decreases and efficiency 
Increases. At a certain pressure the efficiency 
reaches a maximum and then decreases along with the 
specific volume as the pressure continues to Increase. 
This is characteristic of ail recuperated Brayton 
cycles. We will set the pressure at the value corre- 
sponding to the optimum efficiency. In this fashion 
the diameter of the particle and the fractional pres- 
sure drop ore sufficient to specify the specific vol- 
ume and the efficiency. The results are shown in 
Fig. 3(a) for turbine Inlet temperature of 1250° K and 
Fig. 3(b) for turbine Inlet temperature 1150° K. 

The cycle efficiency Is sensitive to turbine Inlet 
temperature. A comparison of Figs. 3(a) and (b) shows 
that the efficiency Increases about four percentage 
points when the turbine Inlet temperature Is Increased 
100 K from 1150 to 1250. 

At a given fractional pressure drop the efficiency 
decreases with Increasing particle diameter because the 
surface area available for heat transfer decreases. 
This, In turn, lowers the recuperator effectiveness 
and, ultimately, the cycle efficiency. 

Generally, the efficiency improves with decreas- 
ing particle size. Also for a given particle size 
efficiency improves with Increasing fractional pres- 
sure drop up to a point. In this range, the recuper- 
ator effectiveness was increasing rapidly enough to 
offset the increased compressor power required to make 
up for that drop. After that point, the recuperator 
effectiveness increase Is not sufficient to overcome 
the Increased losses. 

The fractional pressure drop at this turning 
point is a function of the particle diameter and is 
plotted In Fig. 4 for both 1150 K and 1250 K. The 
optimum value Increases with particle diameter. The 
efficiency value increases with Increasing particle 
diameter, while the specific volume decreases slightly. 
At the larger particle diameters, however, the power 
output Is higher. The fluid velocity, and hence mass 
flow, increases as the size of the particle Increases 


original page is 

POOR QUALITY 


4 


because of the larger entrainment velocity. 


M 

Pr 










For the ahell and tube recuperator with a specl- 
fled fractional pressure drop, tube diameter and flow 
Mach number 'hare la a compressor outlet pressure which 
also maximizes the efficiency. Again we choose that 
pressure as our operating point. The results can be 
seen in Fig. S for a fixed fractional pressure drop of 
.04. The recuperator specific volume Increases with 
Increasing tube diameter for any Mach number and in- 
creases with decreasing Mach number for any tube di- 
ameter. The performance la similar for other pressure 
drops, with the specific volume Increasing with In- 
creasing pressure drop for any tube diameter and Mach 
numbe.. 

On the other hand, efficiency doesn't change with 
Increasing cube diameter at a fixed Mach number and 
pressure drop, it Increases with decreasing Mach number 
for a fixed tube diameter and pressure drop, and It has 
a maximum for a fixed Mach number and tube diameter as 
the pressure drop Increases. When the pressure drop Is 
too low the heat transfer Is too low and not enough 
heat In the turbine exhaust Is recovered. When the 
pressure drop Is too high the compressor work Is too 
large and the Increased heat recovery cannot compen- 
sate. 

Generally, comparing Figs. 3 and 5 one can see 
that for comparable efficiency the fluldlred bed re- 
generators are more compact. This Is the result of the 
more favorable surface to volume ratio of particles 
(6/Dp) compared to the tubes (2/D(). 

CONCLUSIONS 

The high operating temperature capability of flu- 
idized bed regenerators makes them an attractive con- 
cept for advanced cycle application. 


Mach number 
Prandtl number 
f pressure 

Ap pressure drops 

q parameters defined In Eq. (3) and (8) 

H gas constant 

Re Reynolds number 

S surface area defined In Ref. 4 

T temperature 

U heat transfer coefficient In Ref. 4 

V heat exchanger volume 

y ratio of specific heats 

c void fraction 

n efficiency 

U gas viscosity 

C parameters defined In Eq. (7) and (10) 

p density 

T parameters defined in Eq. (7) and (10) 

Subscripts 

c cold particle transfer line 

comp compressor 

g gas 

h hot particle transfer line 

hs heat source 

1 lower bed 

o superficial velocity 








Calculations comparing these fluidized bed regen- 
erators to shell and tube heat exchangers at the same 
conditions with solid density 2.1 grams/cm3 and gas 
velocity of twice the minimum fluidization velocity 
have been made. These calculations show that the total 
volume of the shell and tube heat exchangers is two 
to three times the total volume of the fluidized beds. 
However, the total amount of power which can be han- 
dled by each of these fluidized beds is limited to 
probably less than a megawatt unless very shallow beds 
are considered. The power may be Increased by consid- 
ering rotating beds where the acceleration on the 
solids can be increased above 1 gram (9.8 m/sec2)7,8. 
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APPENDIX - SYMBOLS 
specific heat 

compressor parameter defined In Eq. (3) 
parameters defined In Ref. 4 
diameter 

turbine paramever defined In Eq. 2 
friction factors given In Lq. (22) 
ratio of free board volute to bed volume 
gravity acceleration, 9.806 m^/sec 
mass flow per unit a.ea 
ieat transfer coefflciei t 
bed heights 

length of hot particle transfer line 
mass flow rate 


rec 

recuperator 

8 

solid 

sh 

shell 

t 

tube 

Cl 

terminal 

turb 

turbine 

u 

upper bed 
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igure L. - Schematic diagram of open c^:le gas turbine 
system with fluidized bed heat regenerators. 
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Figure 4. - Fluidized bed recuperator parameters at "optimum" compressor outlet pressure and 
"optimum" fractional pressure drop across the beds. 
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OVERALL THERMODYNAMIC CYCLE EFFICIENCY. V, % 

Figure 5. - The effect of tube diameter and gas flow Mach 
number on the recuperator specific volume and the over- 
all thermodynamic cycle efficiency for shell and tube re- 
cuperator. 
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